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Abstract. 

We compare two fermionic renormalization group methods which have been used 
to investigate the electronic transport properties of one-dimensional metals with two- 
particle interaction (Luttinger liquids) and local inhomogencities. The first one 
is a poor man's method setup to resum "leading-log" divergences of the effective 
transmission at the Fermi momentum. Generically the resulting equations can be 
solved analytically. The second approach is based on the functional renormalization 
group method and leads to a set of differential equations which can only for certain 
setups and in limiting cases be solved analytically, while in general it must be 
integrated numerically. Both methods are claimed to be applicable for inhomogencities 
of arbitrary strength and to capture effects of the two-particle interaction, such as 
interaction dependent exponents, up to leading order. We critically review this for 
the simplest case of a single impurity. While on first glance the poor man's approach 
seems to describe the crossover from the "perfect" to the "open chain fixed point" 
we collect evidence that difficulties may arise close to the "perfect chain fixed point" . 
Due to a subtle relation between the scaling dimensions of the two fixed points this 
becomes apparent only in a detailed analysis. In the functional renormalization group 
method the coupling of the different scattering channels is kept which leads to a better 
description of the underlying physics. 



PACS numbers: 71.10.-w, 71.10.Pm, 73.21.Hb 



Fermionic RG methods for inhomogeneous Luttinger liquids 



2 



1. Introduction 

The sensitivity of one-dimensional (ID) correlated electron systems - so called Luttinger 
liquids (LL) [1] - to local single-particle inhomogeneities (e.g. a single impurity) leading 
to a (back-) scattering of electrons with momentum transfer 2k p (with kp being the 
Fermi momentum) was recognized several decades ago [21 [3l HJ E]. The density- 
density response function at zero frequency shows a power-law singularity when the 
momentum transfer q approaches 2k p with an exponent determined by the Luttinger 
liquid parameter < i^' < 1 (for repulsive interactions; A' = 1 in the noninteracting 
case). It was computed for a specific field-theoretical model - the Tomonaga-Luttinger 
(TL) model - after mapping the interacting fermionic problem onto a free bosonic 
Hamiltonian using the bosonization method [6]. The TL model is known to describe 
the low-energy physics of a large class of ID metals [1]. 

At the beginning of the 90's it became apparent that the realization of sufficiently 
well defined ID quantum wires and thus the possibility to experimentally investigate 
LL physics would only be a matter of a few more years of research. This was one of the 
reasons for a renewed theoretical interest in such systems. The understanding of the 
role of local impurities was further sharpened by renormalization group (RG) studies 
within the bosonized model [3 [8]. In these studies the TL model was supplemented 
in two ways. In a first step a term solely describing the kp to —kp (and vice versa) 
backscattering (of amplitude Vf,) was added to the Hamiltonian. Using bosonization this 
leads to a term oc sin [0(0)], with 0(x) being the bosonic field and a; = the position 
of the impurity. The resulting field-theoretical model is known as the local sine-Gordon 
model. In a RG procedure perturbative in Vf, it was shown that Vf, grows and is thus 
a relevant perturbation (for repulsive interactions with < K < 1). At a sufficiently 
low energy scale the perturbative RG breaks down. In a second step the opposite case 
of two semi-infinite TL chains with a weak link of amplitude twi was studied. This was 
shown to scale to zero. Based on these insights and an exact solution at = 1/2 it was 
speculated that for repulsive interactions even a weak single impurity will scale towards 
the strong impurity limit implying that at low energies observables show power-law 
scaling with exponents characterized by the scaling dimension of the "open chain fixed 
point" 1/K. E.g. the local spectral function p{x ~ 0, u) close to the impurity will vanish 
as uj^^^~^ and the conductance G across the impurity as 



with T being the temperature. For small H the corrections to "perfect" G = e^/h 
conductance are given by the scaling dimension K of the (unstable) "perfect chain 
fixed point", 



Within the local sine-Gordon model the above conjecture was later confirmed using 
numerical [10] and exact [H] methods. 



(1) 



(2) 



Fermionic RG methods for inhomogeneous Luttinger liquids 



3 



The physics of spinless inhomogeneous LLs and systems with spin is quite similar 
in many respects and for simphcity we here mainly focus on the spinless case. The 
scaling exponents just given are slightly different if the spin is included. Towards the 
end of this article we briefly explain how the spin can be included in the two fermionic 
RG approaches described here and discuss a novel feature which results from it. 

In parallel to the bosonization approach a fermionic poor man's RG method was 
developed which does not rely on the mapping onto a bosonic field theory [12]. In this 
approach the single-particle impurity scattering problem is first solved exactly within the 
ID electron gas model. Next, the Hartree and Fock terms [leading order perturbation 
theory in the two-particle interaction U{x)] are computed with the scattering states. 
The Hartree diagram produces a potential term which oscillates as cos {2kpx), decays 
as 1/x, and has a prefactor given by the reflection amplitude of the (bare) impurity. 
Scattering off the Hartree potential leads to a logarithmically diverging correction to 
the transmission amplitude t{k) if k ^ kp 

(3) 



6t{k) ~ —a In 



d{k — kp) 

where ci is a characteristic spatial scale of the interaction and 

„ = M_^ (4) 

is the relevant amplitude of the (Fourier transform of the) interaction, with vp being 
the Fermi velocity. The authors argue that one can iterate this procedure and produce 
divergent terms of higher order [Q;ln{l/|(i(/i; — /ci?)|}]". 

To resum these terms a real space RG method is suggested. At scale A = vpjl 
the scattering potential ~ cos {2kpx)/x is considered only for x G [— /, /] with a spatial 
cutoff /. This way the logarithmic divergence is regularized. Then the region [—1, 1] with 
its renormalized reflection and transmission amplitudes is viewed as an effective barrier 
from which the Hartree potential and the Fock term for the region outside [— /, /] can 
be computed. One next considers Z — > / + 5/ and iterates the above steps. This leads to 
a single differential equation for the transmission amplitude at momentum kp. For the 
following analysis it is essential to note that the nontrivial spatial dependence of the 
scattering potential at each step is ~ cos {2kpx)/x and does not change during the RG 
procedure. 

It was argued that, in contrast to the bosonic RG method which is perturbative 
in the impurity strength, this fermionic procedure, although being applicable only for 
sufficiently small a (weak two-particle interaction), can be used to describe the full 
crossover from the weak to the strong impurity limit. Further down we challenge 
this statement. A detailed analysis for the single impurity case shows that the poor 
man's approach heavily relies on a subtle relation between the scaling dimensions of the 
"perfect" and "open chain fixed points". The method is commonly regarded as being 
very intuitive and has been generalized to study transport through two barriers (resonant 
tunneling) [HKIl], through Y-junctions of three LLs [l5], and through regular networks 
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of LLs [15] as well as transport across a LL-superconductor junction [T7j. Remarkably, 
the lack of a formal framework from which the poor man's flow equations can be deduced 
led to a dispute about their correct form for the double barrier case [I3l [H] . 

In both the fermionic RG as well as the field-theoretical description only the 
scaling of the backscattering amplitude kp — > —kp (and vice versa) was considered. 
Undoubtedly the backscattering component Vf, = Vkp-kp of ^ general scattering 
potential with matrix elements V^^k' is (i) the one which grows the fastest when the 
infrared cutoff A is lowered and is (ii) the component which at the end of the RG flow 
determines the conductance across the impurity. But in more elaborate RG procedures 
all scattering channels — > fc' will couple and also the amplitudes Vk,k' with k, k' not 
equal to -^kp will grow. Even if the matrix elements with momentum transfer different 
from 2kF are absent in the original Hamiltonian generically they will be generated during 
the RG procedure. One might speculate that this feedback will lead to intermediate fixed 
points not captured by the local sine-Gordon model (as no other matrix elements than 
Vkp-kp are kept in the model) and the poor man's RG (as the feedback is neglected). 
Motivated by this reasoning an alternative fermionic RG procedure was developed 
starting from the rather general functional RG (fRG) framework [T9] . It allows for a 
description of the full functional form of the renormalized effective impurity potential. 
Using truncated fRG flow equations it was shown that indeed no intermediate fixed 
points are realized [20l EH [22] . As a further advantage this RG method can directly 
be applied to microscopic lattice models and does not require a mapping onto the ID 
electron gas or an effective field theory. It is rather straightforward to extend the fRG 
approach to study resonant transport through a double barrier [22], transport through 
a wire with three and four impurities [23], and junctions of more than two LLs [241 [25] . 

In the fRG formalism one derives an infinite hierarchy of coupled differential 
flow equations for the cutoff dependent one-particle irreducible vertex functions. In 
practical applications this hierarchy must be truncated. To study the physics of local 
inhomogeneities in spinless LLs up to leading order in the two-particle interaction it is 
sufficient to keep the flow of the frequency independent self-energy (one-particle vertex) 
[20] . The quantitative agreement to the results derived by bosonization in the limits of 
weak and strong impurities as well as to numerical results obtained for small systems can 
be improved if the flow of the static part of the effective two-particle interaction (two- 
particle vertex) is kept [211 122]. Details are given below. The frequency independent 
self-energy can be interpreted as an effective, renormalized impurity potential. For 
strong bare impurities it spatially oscillates with frequency 2k p and decays as the 
inverse distance from the inhomogeneity up to a scale vp/T, beyond which it dies 
off exponentially. Scattering off this effective impurity then leads to the power-law 
decay of the conductance as a function of T. For weak bare impurities and not too 
low temperatures the oscillatory effective potential only decays with exponent K. For 
lower temperatures a crossover to the strong impurity behavior can be observed and the 
system flows to the "open chain fixed point" . 

Below we give more details on the truncated fRG procedure and the basic equations 
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for the poor man's RG (Sec. [2]). We discuss the most important results and compare the 
outcomes of the two fermionic RG approaches (Sec. [3]). It is shown that the truncated 
fRG can be used to describe the full crossover from the weak to the strong impurity 
limits while the poor man's RG seems to have difficulties close to the "perfect chain 
fixed point". We show that a single RG equation - as derived in the poor man's RG 
for the effective transmission at momentum fc^;' - is in general not enough to describe 
the crossover behavior. We then briefly comment on how the spin degree of freedom 
is included in both methods (Sec. Hj). Finally, the comparison of the two fermionic RG 
methods for the single impurity case is summed up and briefly extended to Y-j unctions of 
three LLs (Sec. [5]). In the Appendix we discuss the additional approximations necessary 
to derive the poor man's RG equation within the fRG framework. 

2. The methods 

In this section we present the basic equations of the fRG approach and the poor man's 
RG. We focus on the simplest situation of a single impurity in an otherwise "perfect" 
chain which is coupled "adiabatically" to two semi-inflnite noninteracting leads [9]. 

2.1. The functional RG approach 

The fRG was recently introduced as a powerful new tool for studying interacting Fermi 
systems [18]. It provides a systematic way of resumming competing instabilities and 
goes beyond simple perturbation theory even in problems which are not plagued by 
infrared divergences. In our application the interacting quantum wire will be coupled to 
noninteracting leads as it is the case in systems which can be realized in experiments. 
The leads are modeled as ID tight-binding chains with hopping matrix element t = 1. 
As we are interested in the effect of a single impurity in the bulk part of the wire and not 
the role of the contacts the latter are modeled as being perfect, that is the conductance 
through the interacting system of length in the absence of the impurity is e'^/h (up to 
our numerical accuracy). This implies that the interaction must be turned on spatially 
smoothly close to the contacts ( "adiabatic contacts" ) [221 12] • Before setting up the fRG 
scheme we integrate out the noninteracting leads using a standard projection technique 
so that we only have to deal with an interacting system of size A^ [22] . The energy scale 
6 = vp/N related to the length of the interacting part will provide an infrared cutoff 
for power-law scaling as a function of other energy variables. 

In the fRG procedure the noninteracting propagator Qq is replaced by a propagator 
depending on an infrared cutoff A. Speciflcally, we use 

QQ{iUn) = X^{^^n)Go{iuJn) (5) 

with A running from oo down to 0. The cutoff function is chosen such that the 
fermionic Matsubara frequencies |a;„| ^ A are cut out (for details see the Appendix A 
of Ref. [2^). Note that Qq contains the frequency dependent self-energy contribution 
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from the leads 



= !^ ^1 _ ^1 _ j S,,,, + , (6) 

with G [1,A^] indicating the lattice site and the filling n-, interaction U-, and T- 
dependent chemical potential /z = fi{n,U,T). It is chosen such that the entire system 
(wire and leads) has the desired filling at any U and T. Using in the generating 
functional of the irreducible vertex functions and taking the derivative with respect to 
A one can derive an exact, infinite hierarchy of coupled differential equations for the 
vertex functions, such as the self-energy and the irreducible two-particle interaction 
r^. In particular, the fiow of the self-energy (one-particle vertex) is determined by T,^ 
itself and the two-particle vertex, while the fiow of is determined by S^, F^, and the 
fiowing three-particle vertex Fg . The latter is computed from a fiow equation involving 
the four-particle vertex, and so on. At the end of the fRG fiow is the self-energy 

S of the original, cutoff-free problem we are interested in. A detailed derivation of the 
fRG flow equations for a general quantum many-body problem which only requires a 
basic knowledge of the functional integral approach to many-particle physics and the 
application of the method for a simple problem are presented in Ref. [T9] . 

In practical applications the hierarchy of flow equations has to be truncated and 
^A=o Qj^iy provides an approximation for the exact S. As a flrst approximation we 
neglect the three-particle vertex. The contribution of Fg to F^ is small as long as F'^ 
is small, because Fg is initially (at A = oo) zero and is generated only from terms of 
third order in F^. Furthermore, F^ stays small for all A if the bare interaction is not 
too large. This approximation leads to a closed set of equations for F^ and [22j. 

We here do not give these equations but instead implement a second approximation: 
the frequency dependent flow of the renormalized two-particle vertex F^ is replaced 
by its value at vanishing external frequencies, such that F^ and hence remains 
frequency independent. Since the bare interaction is frequency independent, neglecting 
the frequency dependence leads to errors only at second order for the self-energy, and 
at third order for the two-particle vertex at zero frequency. For the approximate flow 
equations we then obtain 

S^, , = E E ^^^'(^^n) 2,1,2 (7) 



dA 2-K 

|a;„|ssA 2,2' 



and 



^r^',2';l,2 — ^ E/ E/ I 2 ^sls'l^^") ^^n)ri'',2';3,4 r 



A 

3',4';1,2 



|a;„|RiA 3,3',4,4' 



+ ^3^3'(^^n) ^M'(^^n) [ ^f/ 4/.]^,g Fg, 2/.4_2 + r2',4';l,3 r3/,l';4,2] | ) (8) 

where the lower indexes 1,2, etc. stand for the single-particle quantum numbers of the 
chosen basis (e.g. the Wannier basis or the momentum state basis) and 

G\lUn)=[Q,\li^n)-T.''y' . (9) 
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The abbreviation \ujn\ ~ A stands for taking the positive as well as negative Matsubara 
frequency with absolute value closest to A at fixed T. At the initial cutoff A = oo the 
flowing two-particle vertex Ty 2'-i 2 given by the antisymmetrized interaction and the 
self-energy by the single-particle terms of the Hamiltonian not included in Qq, e.g. the 
impurity. 

Next we specify our model for the interacting LL wire with noninteracting leads. 
We use the ID tight-binding model with nearest-neighbor hopping t = 1 and nearest- 
neighbor interaction Ujj+i. It is given by 

H = ifkin + -f^int 

00 N-1 

j=-oo j=l 

where we used standard second-quantized notation with c] and c- being creation and 
annihilation operators on site j, respectively and the local density operator uj = Cj Cj. 
The filling n and interaction dependent real number p must be determined iteratively 
such that at fixed common n (for the leads and the wire) the interacting wire has the 
desired filling. For the mostly studied case of half-filling particle-hole symmetry implies 
u = 1/2. As mentioned above the interaction f/jj+i is turned on and off smoothly over 
a small fraction of the interacting lattice sites. The constant bulk value is denoted 
by U. 

The homogeneous model with a constant interaction U across all bonds can be 
solved exactly by the Bethe ansatz p7] . It shows LL behavior for all particle densities n 
and any interaction strength except at half filling for \U\ > 2. The U- and n-dependent 
LL parameter K can be determined solving coupled integral equations j2Tj, which in 
the half-filled case can be done analytically with the result 

K = \- aiccos ( -^\\ =1-- + 0{U^), (11) 
vr V 2 / J vr 

for \U\ < 2. As depends on the filling it is essential to tune the interacting wire to 
the desired filling by choosing the appropriate u [see Eq. (ITOl) ]. Only then a comparison 
of the fRG results to the ones obtained by bosonization is possible. 
The Hamiltonian ffTOl) is supplemented by an impurity part 

^i-p = E^..'4s" (12) 

where Vjj' is a static potential. Local site impurities of amplitude V at jo are given by 
a potential 

^jd' = ^^jd'^j,jo (13) 
and local hopping impurities with hopping t' across the link jo, jo + 1 by 

^,/ = ^',. = -(^'-l)Wi^.,.o- (14) 
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In our RG procedure we parameterize the static real space vertex by a flowing 
nearest-neighbor interaction 

We furthermore neglect the feedback of the self-energy on the flow of U^. The differential 
equation for the renormalized bulk interaction is then determined after transforming 
Eq. (fTSl) to momentum space and fixing the momenta on the Fermi surface [28]. It reads 
(in the limit —>■ oo) 

= h{Co^) {U^f (16) 

with 

1 vi 

h{x) = — Re -{fi + ix 

2n 




"i,! -I- 1 Rtii -I- (Mt"^ -I- i-r'^ l 

(17) 



3zyU^ — IQjTx — 12ifi^{x^ + 1) + Qx^fJ- + 18x/i + Qix^ + ix 

X 



7r(2/i + ix){A — fi'^ + x'^ — 2ixjiY 
where cj„ stands for the fermionic Matsubara frequency closest to A. The initial 
condition is limAo^oo = U. The bond dependent renormalized interaction f^/j-ti 
follows from the bulk by multiplying with Ujj±i/U. After performing these 
additional approximations the flow equations for the Wannier basis matrix elements 
of the self-energy read 

^ ^i^i = ~ 27r ^ ^ ^^^^^ ^t+rj+ri^^-^) ' 
|<^„|«A)-=±1 

^^h^i = ut,^i0tm(^^n) . (18) 

|(^„|»A 

All other components vanish. The self-energy enters also the right-hand side of these 
equations via the full propagator Q^. The initial conditions at A = Aq ^ cxd are 

S^Ao = [1/2 - i^in, U)] (f/,_i,, + f/,,,+i) , S^Ao^^ = . (19) 

At the end of the fRG procedure the diagonal part of the approximate self-energy j can 
be interpreted as an effective impurity potential, while — 1 + ^j,j±i is the renormalized 
hopping. 

For general parameters the set of coupled differential equations can only be solved 
numerically. This can be done for fairly long chains (up to 10^ lattice sites) as computing 
the right-hand side of the flow equations only requires the tridiagonal part of the inverse 
of a tridiagonal N x N matrix which can be obtained in order time [28]. In the 
limiting cases of weak and strong impurities and neglecting the flow of the two-particle 
interaction altogether the flow equations can be solved analytically at T = 0. This is 
discussed in Sec. [31 

The truncated fRG procedure provides an approximation for the self-energy and 
thus via the Dyson equation an approximation for the interacting one-particle Green 
function Q. If the inelastic (frequency dependent) processes the two-particle interaction 



Fermionic RG methods for inhomogeneous Luttinger liquids 



9 



generates are neglected, the conductance can be computed from the (1, A^) matrix 
element of the retarded Green function Q"^ by a generalized Landauer-Biittiker formula 

m 

G{T,6) = jf^'^ (-£) \tis,T,6)\'de, (20) 

where f{e) = l/(e^^~^^/^ + l) is the Fermi function, B the band width (5 = 4 for t = 1), 
and 

|t(£,,T,5)|2 = 4sin2(fc) \gij,{ek)\' (21) 

the effective transmission. Because of the frequency independence of S the analytic 
continuation of Q to the real axis, necessary to obtain Q^{e), is straightforward. 
During the RG procedure the self-energy and thus the propagator acquire additional 
dependencies on the temperature and on the energy scale 6 ^ 1/N. 

Inelastic processes vanish in the limit T ^ 0. For T > they are absent in our 
approach as we assumed the self-energy to be frequency independent and using Eq. fl20l) 
is thus consistent with our earlier approximations - it presents no further approximation. 
In general inelastic processes lead to corrections to Eq. (120|) of order U"^ and higher. It 
is important to note that inelastic processes due to the two-particle interaction are also 
neglected in the poor man's RG as well as in the bosonization approach. In the latter 
they are disregarded because the bulk part of the local sine-Gordon model corresponds 
to the TL model in which all infrared RG irrelevant terms such as umklapp scattering 
are left out. At finite scales such terms imply inelastic scattering. 



2.2. Solution of the poor man's RG equation 

The basic idea of the poor man's RG was already presented in the introduction. It is 
set up for the ID electron gas model. The method leads to the flow equation [12] 

"^^"^ at^(l-k^l') (22) 



d\n{A/[vF/d]) 

for the transmission t^ at kp, with d and a introduced in Sec. [H It only depends on 
a single energy scale, the infrared energy cutoff A, which must be contrasted to the 
fRG approach in which the scales A, T, and 6 appear. The integration of Eq. (l22ll 
starts at the ultraviolet cutoff Aq = vp/d with the initial condition t^" = t^, where 
to is the bare transmission. To obtain the energy or temperature dependence of the 
renormalized transmission the integration is stopped at A = e or T, respectively. The 
assumed dependence on only one energy scale implies that the variable one is interested 
in must be much larger than all the other scales. In the fRG in contrast the interplay 
of the scales A, T, and 6 can be studied. The solution of Eq. fl22|) is 
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where s must be substituted by the energy variable of interest (e.g. s T). The bare 
reflection probability is denoted by |roP = 1 — |toP- The expression for the conductance 
follows after multiplying the transmission (probability) by e'^/h. 

For the upcoming analysis it is important to note that the leading order contribution 
of K for the ID electron gas model with two-particle interaction U{x) is given by [1] 

U{0) U{2kF) 



K 



2nvF 2ttvf 



a 



(24) 



with the forward- [?7(0)] and backward-scattering [C/(2fci?)] amplitudes. 

In the Appendix we derive the poor man's RG equation ( |22l) for the ID electron 
gas model within the fRG framework. This requires approximations which go beyond 
the ones used to obtain Eqs. ([7]) and (IHl). 



3. A comparison of results for the conductance and the renormalized 
impurity potential 

In this section we compare the results for the conductance and the renormalized impurity 
potential obtained using the two fermionic RG methods. We start with the two limits 
of weak and strong bare impurities. 



3.1. The conductance for weak and strong impurities 



The solution Eq. 
limits of weak (|ro/t, 
this gives 



of the poor man's RG equation can easily be expanded in the 
<^ 1) and strong (|to/''^o| ^ 1) bare impurities. In the first case 



.s|2 



1 





2 


S 


to 




Ao 



-2a 



while in the second one finds 



J.SI2 



to 


2 


S 


^0 




To 



2a 



(25) 



(26) 



Within the bosonization approach the two relevant scaling exponents are give by 2{K—1) 
and 2(1/ K - 1) [see Eqs. and With Eq. ^ one finds 

2{l/K-l)^-2{K-l)^2a (27) 

to leading order in the two-particle interaction and Eqs. (125|) and fl26|) are consistent 
with the bosonization results Eqs. ([2]) and ([T]) if s is replaced by T. One is thus tempted 
to conclude that the poor man's RG result Eq. fl23p describes the full crossover from the 
"perfect chain fixed point" to the "open chain fixed point" as long as the interaction 
remains small a <C 1. Below we argue that the analytic Breit-Wigner form Eq. fl2^ is 
incomplete and that the apparent agreement with the bosonization result heavily relies 
on the subtle relation Eq. (1271) (which can be traced back to what is called "duality" in 
the field-theoretical approach). 
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In the two limits considered in the present subsection also the fRG flow equations 
can be tackled analytically provided that the flow of the effective nearest-neighbor 
interaction is neglected. In Ref. [20] it was shown that at T = for weak bare impurities 

„fc ~ A-[U{0)-U{2kF)]/{2nVF) ^28) 

rather independently of the model considered (that is not only for the lattice model 
with nearest-neighbor interaction and independent of the type of impurity considered). 
Within the Born approximation of single-particle scattering theory the transmission 
(conductance) can be obtained from the backscattering amplitude S^^ Assuming 
that A in Eq. (128|) can be replaced by the scale s (as it is routinely done in the poor man's 
RG; see above) the transmission shows the correct scaling behavior with an exponent 
which agrees to the bosonization result to leading order in the two-particle interaction. 
In Ref. [20] it was also shown that at T = 0, for the lattice model described above, 
and for a strong hopping impurity (weak link) the local spectral function close to the 
impurity shows power-law scaling (as a function of energy) with an exponent equal to 
that found for an open chain. The latter agrees to leading order with the bosonization 
exponent 1/K — 1. Following standard arguments (Fermi's Golden Rule) [7] one can 
then conclude that the conductance at low energy scales shows power-law scaling with 
an exponent which agrees to leading order with the bosonization result 2(1/ K — 1). 

3.2. The generalized Breit-Wigner form of the transmission 

The expressions ( l20l) and (|2T1) for the transmission (conductance) used at the end of the 
fRG procedure can be brought into a form similar to that of the poor man's RG Eq. fl23|) . 
To be speciflc we consider a site impurity at Jq (in the bulk of the interacting wire) of 
amplitude V. In a flrst step one runs the fRG flow down to A = and determines the 
approximate self-energy S. Denoting the single-particle version of the noninteracting 
part of the Hamiltonian by ho, the N x N matrix 

h = ho + J: + (29) 

forms an effective single-particle "Hamiltonian" , which is non-hermitian because of the 
leads contribution S'^^^ Eq. (I6|) . After setting the renormalized hopping matrix element 
t = — 1 + from site jo — 1 to site jo to zero, that is after cutting the chain in two 

parts, we deflne an auxiliary Green function [resolvent matrix of size {N —jo) x (N — jo)] 

g{z) = [z-h{i=0)]'^ . (30) 

Via S the auxiliary Green function still contains information about all system 
parameters, in particular T and 6: Q{z,T,6) = G{z). Within the framework of single- 
particle scattering theory and using a standard projection technique [22] the effective 
transmission can be written as 

\tie,T,6)\' = l£!(£ili^) ^ (31) 

e -V{T,5) -2n{e,T,6) +AT^{e,T,6) 



Fermionic RG methods for inhomogeneous Luttinger liquids 



12 



with 

POj,-i,j,-iie + tO, T, 5) = n{e, T, 5) - iT{e, T, 5) , (32) 

real functions Q and F, and the renormahzed energy of the impurity site jo, V{T, 6) = 
V + 

^iojo(^''^)- ■'■f '^'^^ consider T — and plug Eq. ([31]) into the generalized 
Landauer-Biittiker formula fl20|) we end up with 



4gH<(0.0..)|-, . (33) 

e 



V{0,6) + 2n{0,0,6) +4F2(0,0,5) 

Using s = 5 as our scaling variable this generalized Breit-Wigner form is similar to 
that of Eq. fl23|) obtained within the poor man's RG. Using the energy 6 related to the 
inverse size of the interacting part of the system as an energy scale is valid also in the 
poor man's RG according to the construction of the RG equation. The characteristic 
difference of the two approaches is the absence of a scale dependence in the first term 
of the denominator of Eq. (I23l) . This will become crucial in the following. 

To further develop the relation between the two generalized Breit-Wigner forms we 
next analyze the scaling properties of Q and F. In Fig.[T]results for V{S)+2Q{6) and F(5) 
obtained by numerically solving the truncated fRG flow equations for strong [Fig. [T](a)] 
and weak [Fig. [It^b)] bare impurities are shown on a log-log scale. The parameters are 
U = 0.5 and 1, n = 1/2, 1/ = 20 [Fig. W^a)] and V = 2 ■ IQ-^ [Fig. [I](b)], and the size 
of the interacting wire varies from = 257 to = 65537. For strong bare impurities 
V{S) + 2Q{S) is constant, while the imaginary part of the auxiliary Green function shows 
power-law scaling as a function of 6. As T{S)/[V{S) + 2Q{6)] <^ 1 the imaginary part 
can be neglected in the denominator of Eq. (1331) . From this one concludes that the 
exponent of G{6) is twice as large as that of F(5). The exponent extracted from the 
fRG data for T{6) [see Fig. [11(a)] agrees to leading order in U with the expected result 
1/K — 1, with K taken from Eq. (fTTl) [22] . For weak bare impurities T{6) is constant 
while V{6) + 2n{S) follows a power law. Furthermore, [V{6) + 2n{6)]/T{6) < 1 and 
the conductance Eq. (l33l) can be expanded in this ratio. From this expansion one 
concludes that the correction of the conductance to e'^/h scales with twice the exponent 
of V{6) + 2Q{6). The exponent of V{6) + 2Q{6) [see Fig. [H^b)] determined numerically 
indeed agrees to leading order with the expected one K — 1 [22] . 

It is important to note that for T = Eq. (l33ll provides the exact analytical form of 
the conductance even if one goes beyond the fRG truncation used here. We expect that 
if computed exactly V{6) + 2Q{6) and T{6) show power-law scaling with the bosonization 
exponents (not only the leading order exponents) in the corresponding limits of weak 
and strong impurities. 

Identifying 4T^{0,0,S) ^ |toP|5/AoP° and [V{6) + 2^(5)]^ ^ |rop might appear 
ambiguous as one can devide the numerator and denominator of Eq. f l33|) by the scale 
dependent part of [V^(5)+2fi(5)]2 and end up with a function similar to Eq. f[23l) in which 
the first term in the denominator is scale independent. However, the above relations 
become unique if one considers the finite e generalization of Eq. ([23]) as given in Eq. (5) 
of Ref. [I3j which must be compared to Eq. ( [311) . 
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Figure 1. 6 dependence oi V + (solid blue) and F (dashed red) entering the 
generalized Brcit-Wigner form Eq. ([55]) aX n — 1/2, U = 0.5 (circles) and [7 = 1 
(squares) obtained from the truncated £RG. In (a) the case of a strong site impurity 
with y = 20 is shown, while in (b) a weak impurity V — 2 ■ 10^^ is considered. 



These insights provide hints that the absence of the scale dependence in the first 
term of the denominator of Eq. (l23l) (the "real part") constitutes a problem for the 
poor man's RG in the limit of weak bare impurities. In fact, in the poor man's RG 
approach the "imaginary part" carries the scaling in both limits of weak and strong 
impurities. The observation that Eq. fl25]) reproduces the bosonization results in both 
limits thus relies heavily on the fact that 1/ K — 1 = —{K — 1) to leading order in the 
two-particle interaction. This is different in the fRG approach, in which two different 
functions carry the scaling properties in the two limits. The analysis also reveals that 
for a method that leads to a generalized Breit-Wigner form one needs at least two 
independent RG equations (one for the real and one for the imaginary part of ^) to 
fully describe the crossover. 

3.3. The one-parameter scaling function 

We now further elaborate on the differences of the two fermionic RG approaches by 
studying the one-parameter scaling properties of the conductance (transmission). 

That the full crossover from the weak to the strong impurity limit within the local 
sine-Gordon model follows a one-parameter scaling function was shown analytically in 
Ref. for = 1/2. Within this model it was later confirmed for other LL parameters 
plot [TT] and generalized to microscopic models in Refs. [21] and [22] . The scaling ansatz 
is given by [7] 

G{s) = ^Gk{x) , X = [s/soiU, n, V)f-' (34) 
n 

with the nonuniversal scale sq and s being T or 6. For appropriately chosen sq the 
G{s) curves for different V (but fixed K) can be collapsed onto the i^-dependent 
scaling function Gk{x). It has the limiting behavior Gk{x) ~ 1 — for x ^ and 
Gk{x) ~ x'"^/^ for X oo. For the following it is important to note that using either 




Figure 2. One-parameter scaling plot of the conductance obtained by the truncated 
fRG procedure. Open symbols represent results obtained for U = 0.5, n = 1/2, and 
different T and V, while filled symbols were calculated for U = 0.851, n — 1/4. Both 
pairs of U and n lead to the same K^'^ — 1 = —0.15. The solid red lines indicate the 
asymptotic behavior for small and large x. The solid black line indicates a power law 
with exponent —2. 



Eq. ffTTj) or Eq. fl24l) the large x exponent —2/K has a contribution that is hnear in the 
two-particle interaction. It follows from taking the "open chain fixed point" exponent 
2(1/ K — 1) (large x behavior) and dividing by half the exponent K — 1 of the "perfect 
chain fixed point" which by convention was explicitly taken into account in the definition 
of the variable x [see Eq. flM|) ]. 

Achieving the form Eq. ( !34l) for the result of the poor man's RG Eq. ( !23l) is 
straightforward. The scaling variable is x = \ro/to\{s/Xo)^°' and the scaling function 

GTi^) = r:^ • (35) 

Independent of the strength of the two-particle interaction it is thus given by the 
noninteracting scaling function Gk=i- While the behavior at small x is reproduced 
correctly, for large x, G^{x) decays ~ This can be traced back to the fact 

that within the poor man's RG the weak and strong impurity exponents are both 
strictly linear in the two-particle interaction and have the same modulus. As mentioned 
above the correct large x exponent —2/K has a contribution linear in the two-particle 
interaction. 

In Fig. [2] we show conductance data obtained numerically from the fRG procedure 
at a fixed small 6 and for temperatures T ^ 6. For appropriately chosen so{U, n, V) the 
data for different V and U = 0.5, n = 1/2 (open symbols) fall onto a single curve. The 
different symbols and colors indicate different V. A similar collapse is found for T = 
and s — *• 5 [21]. In the definition of x [see Eq. (IMll ] we used the fRG approximation for 
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K -1 {ioiU = 0.5 and n = 1/2, K^'^ - 1 = -0.15). The same - 1 can be found at 
filling n = 1/A for U = 0.85. Data for these parameters are shown as the filled symbols 
in Fig. [21 The collapse confirms that the one-parameter scaling function indeed only 
depends on K and not on U and n separately. 

In contrast to G^{x), G^'^{x) depends on K and for large x shows a power- law 
decay with an exponent different from —2 (see the black line in Fig. [2]). In the fRG 
approximation the two exponents 1/K—l and —{K—1) are determined numerically from 
the temperature dependence of the conductance for weak and strong bare impurities (or, 
with the same result, as described in the last subsection). In contrast to the poor man's 
method they have higher (than leading) order corrections and the [/-dependent part 
does not cancel if the scaling ansatz is applied. In fact, as shown in Figs. 5 and 7 of 
Ref. [22] the bosonization and fRG exponents agree well beyond the linear regime in 
the limits of weak and strong impurities, even though formally only the leading order 
Taylor coefficients are the same (as terms of order U"^ and higher are only partially kept 
in the truncated fRG scheme). In the light of these observations it is interesting to 
ask whether the fRG approximation of the ratio 2(1//^ — 1)/{K — 1) = —2/K, that is 
the large x decay exponent of the scaling function, has the correct leading order (in U) 
behavior. Unfortunately, the numerical accuracy of the data for the two exponents is 
not high enough to make definite statements about this. For small U (that is in the 
linear regime of —2/K), both exponents are small and any small error (in particular 
of the denominator) will lead to an inaccurate ratio. For K close to 1, for which the 
truncated fRG works best, the complete form of the scaling function (beyond the two 
limits X and x oo) is not known analytically within the local sine-Gordon 
model, and no direct comparison with the fRG approximation is possible. But even for 
K = 1/2, the largest LL parameter (that is the smallest two-particle interaction) for 
which Gk{x) was computed in the field-theoretical model, the fRG data of our lattice 
model are surprisingly close to the curve obtained using bosonization as shown in Fig. 4 
of Ref. [21]. In this respect the fRG method goes significantly beyond the poor man's 
approach. 

3.4- The decay of the effective, renormalized impurity potential 

It seems to be possible to understand the problems of the poor man's RG close to the 
"perfect chain fixed point" in terms of the real space properties of the effective impurity 
potential. As described in the introduction in the construction of the poor man's RG 
equation the effective impurity potential always decays as the inverse distance from the 
impurity, while the prefactor is renormalized. This holds up to the spatial scale set 
by the infrared cutoff (e.g. the temperature T) and regardless of the strength of the 
bare impurity. Using the fRG approach we now show that this is correct for strong 
bare impurities and on asymptotically small scales. For weak bare impurities and on 
large to intermediate scales the renormalized impurity potential decays more slowly, an 
observation which can also be related to the results known from bosonization. 
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Figure 3. Decay of the oscillatory part of the off-diagonal matrix element of the self- 
energy away from a single hopping impurity at bond jo, jo + 1. Results for t' = 0.1, 
jo = 5000, N = 10"^, U =l,n= 1/2 and different temperatures T = 10"^ (sohd red), 
T = 10-2 (short dashed blue), T = lO^^ (dashed-dotted green), and T = IQ-^ (long 
dashed magenta) are presented. The left panel shows the data on a log-log scale, the 
right panel on a linear-log scale. For comparison the left panel contains a power law 
(j ~ jo)~^ (thin black line). 

In Fig. [3] fRG results for the decay of the oscillatory part of ^jj+i away from jo 
for a strong bare impurity at different T are shown. Here we used a hopping impurity 
with t' = 0.1, that is a hopping across the bond jo, jo + 1 which is only 10% of the 
hoppings across all other bonds. Generically, for the lattice model under consideration, 
both the diagonal and off-diagonal parts of E show decaying oscillations. For symmetry 
reasons (particle- hole symmetry) Sjj = for half- filling n = 1/2 and hopping impurities, 
as chosen in Fig. [31 The figure clearly shows that ^jj+i ~ |j — jo|^^ up to a scale 
~ T^^ beyond which the oscillation decays exponentially. For more general parameters 
and strong bare impurities Sjj ^ shows the same characteristics. This result is in 
agreement with the reasoning leading to the poor man's RG equation. 

In Fig. m we show the effective exponent for the decay of the oscillations of E^j 
as a function of the distance from a site impurity of different strengths V obtained at 
T = 10~^ (dashed) and T = 10~^ (solid). For weak bare impurities and on intermediate 
scales the effective impurity potential follows a power- law decay with a fZ-dependent 
effective exponent (see the plateau at value ~ 0.71 for V = 0.01 at T = 10"''). 
Only for larger distances and lower temperatures, that is on lower-energy scales, or 
for stronger bare impurities this behavior crosses over to the strong impurity limit with 
decay exponent 1. This behavior seems not to be included in the poor man's approach. 
The steep rise of the effective exponent on a scale j — jo ^ Vp/T indicates the crossover 
to exponential decay (see Fig. [3]). The off-diagonal component of the self-energy 
shows similar behavior. 

When studied as a function of U the effective exponent at the weak impurity 
plateau is found to be close to K. This exponent K can be related to the results 
known from bosonization. For weak impurities and on large to intermediate energy 
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Figure 4. Effective exponent for tiie decay of oscillations of 'Ejj as a function of the 
distance from a site impurity of strengths V = 0.01, 0.1, 0.3, 1, 3, 10 (from bottom to 
top) at r = 10^3 (dashed) and T = 10"^ (sohd), half-fiUing n = 1/2 and interaction 
strength U = 1. The impurity is located at the center of a chain with N = 10^ sites. 

scales the effective impurity stays small and is approximately ~ Q{vf/T — \j — 
jo\) cos[2kp{j — jQ)]/\j — jo\^ ■ To obtain the reduction of the transmission with 
respect to 1, one can then treat the renormalized impurity potential within the Born 
approximation. This requires the 2kp component of the impurity potential. Fourier 
transforming the above potential leads to S ^ and finally to G 6 //i ~ 

_j'2{K-i)^ that is the scaling behavior known from the local sine-Gordon model. 

4. Including the spin degree of freedom 

Both fermionic RG methods were extended to study inhomogeneous LLs with spin. 
Within the generalization of the local sine-Gordon model to the case with spin the 
physics is in many respects similar to that of spinless inhomogeneous LL. The system is 
characterized by the "perfect" and the "open chain fixed points" with scaling dimensions 
K/2 + 1/2 and 1/{2K) -1-1/2, respectively [3, [H]. We here assume a spin-rotationally 
invariant interaction. From the scaling dimensions the exponents of the power-law 
corrections to the fixed-point conductance and the exponent of the local spectral weight 
follow as in the spinless case. 

To obtain the scaling exponents to first order in the two-particle interaction within 
fermionic RG methods it is essential to keep the flow of the static two-particle vertex (the 
effective interaction). As it is known from the so called g-ology analysis for homogeneous 
LLs the backscattering component between electrons with opposite spin flows to zero 
if the infrared cutoff is lowered [29j. In the poor man's RG the flow equation for the 
transmission at momentum kp was supplemented by the g-ology RG equations for the 
static two-particle vertex of a homogeneous TL model. This way the prefactor a in 
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Figure 5. Temperature dependence of the conductance for the extended Hubbard 
model with = 10'' sites and a single site impurity of strength = 10, for a Hubbard 
interaction U — I and various choices of the nearest-neighbor interaction U' obtained 
by fRG. The density is n — 1/2, except for the lowest curve, which has been obtained 
for n = 3/4 and U' = 0.65 (leading to a very small backscattering interaction). The 
dashed line is a power-law fit for the latter parameter set. 



Eq. (|22l) becomes A- dependent. For strong bare impurities a new feature appears. If for 
the two-particle interaction U{0) — 2U{2kF) < 0, that is for repulsive interactions with 
a sizable bare backscattering component, the conductance across a single impurity first 
increases as a function of the energy scale s before it is suppressed on asymptotically 
small scales. This effect is not captured within the local sine-Gordon model. 

The fRG approach was generalized to the extended Hubbard model with local 
interaction U and nearest-neighbor interaction U' [26]. In a wide range of parameters 
- to which we stick in the present context - this model shows LL behavior. The 
flow equations for the different components of the static two-particle vertex go beyond 
the g-ology approximation (see Figs. 2 and 3 of Ref. [26]), but the feedback of the 
inhomogeneity on the flow of the vertex (via the self-energy) is neglected. As shown 
in Fig. [5] also within the fRG approach for a strong site impurity \^ = 10 the 
conductance first increases when the temperature is lowered before the asymptotic 
power-law suppression sets in. For the extended Hubbard model, U{0) — 2[/(2/cf) = 
2f/'[l — 2 cos{2kp)] — U, which can be positive or negative for U, U' > 0, depending on the 
density and the relative strength of the two interaction parameters. At quarter-filling 
n = 1/2, f/(0) — 2U{2kp) is negative and therefore leads to an enhanced conductance 
for U' < U/2. In case of a negative U{0) — 2U{2kF) the crossover scale is given by 

Tn oc exp I — ^ I , (36) 

\U{{))-2U{2kF)J 

and becomes exponentially small for weak interactions. For f/' = the conductance 
increases as a function of decreasing T down to the lowest temperatures shown in Fig. O 
For increasing nearest-neighbor interactions U' a suppression of G{T) at low T becomes 
visible, but in all the data obtained at quarter- filling n = 1/2 it is much less pronounced 
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Figure 6. Spectral weight at tiie Fermi level near a hopping impurity t' = 0.5 as a 
function of S = vp/N for the extended Hubbard model with U' = U/\/2, for various 
choices oi U. (a,) n — 3/A (leading to small backscattering). (b) n — 1/2 (leading to 
sizable backscattering) . Results from the fRG (open symbols) are compared to DMRG 
data (filled symbols). 



than what one expects from the asymptotic power law with exponent 1/K — 1. By 
contrast, the suppression is much stronger and follows the expected power law more 
closely if parameters are chosen such that the bare two-particle backscattering becomes 
negligible, as can be seen from the conductance curve for n = 3/4 and U' = 0.65 
(magenta) in Fig. O The value of K [SB] for these parameters almost coincides with the 
one for another parameter set in the plot, n = 1/2 and U' = 0.75 (blue triangles), but 
the behavior of G{T) is completely different. Note that at T ~ 5 the scaling is cut off as 
can be seen at the low T end of the magenta curve in the figure. The scaling exponents 
(if accessible, as e.g. for the magenta curve) extracted from the T- or 5-dependence of 
the conductance (for weak and strong bare impurities) as well as the exponent extracted 
from the w-dependence of the local spectral function agree to leading order in the two- 
particle interaction with the bosonization results. 

Using the fRG we are in a position to even address very subtle questions such as 
how accurately the crossover scale is reproduced by the truncated procedure. To this 
end we need exact or numerically exact data to compare our fRG results with. As no 
such data are available for the conductance of the extended Hubbard model with a local 
impurity we compute the spectral weight w at T = 0, close to the impurity, and at 
c<j = /i as a function of 6 for this model and compare it to very accurate numerically 
obtained density-matrix renormalization group (DMRG) data [30]. We here decouple 
the leads (as they are not included in the standard DMRG approach) and study a 
finite system of length A^. In Fig. [6] we present fRG and DMRG results for w taken 
at site Jo for a hopping impurity with t' = 0.5, U' = different U, and n = 3/4 

[Fig. M^a); small backscattering] as well as n = 1/2 [Fig. [6](b); sizable backscattering]. 
The spectral function Pj^u) close to the impurity site and the spectral weight w{6) show 
a similar crossover behavior as G{T) [see Fig. [6](b)]. For parameters at which the bare 
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two-particle backscattering is small as in Fig. E](a) the fRG and DMRG data agree quite 
well up to the smallest S, that is the largest system sizes N that can be studied by 
DMRG. A similar agreement is found in the spinless case. For 5 — 0, w{S) follows a 
power law with an exponent which agrees to leading order in the two-particle interaction 
with the bosonization result 1/{2K) — 1/2. In contrast, for parameters with a sizable 
backscattering as in Fig. [6](b) the DMRG and fRG data agree only for sufficiently large 
6 (small A^). For larger U and U' the deviation increases and the crossover scale of 
the truncated fRG is easily an order of magnitude too small [see the data for [/ = 2 
in Fig. [6](b)]. As for the conductance the power-law decay sets in only for energies 
sufficiently smaller than the crossover scale. We thus conclude that for a quantitative 
description of the crossover behavior specific to the spinful case one further has to 
improve the fRG procedure. We expect that including the feedback of the self-energy 
on the flow of the two-particle vertex, that is the feedback of the impurity on the flow 
of the effective two-particle interaction, is essential here and will lead to an improved 
approximation. The fRG provides a formal framework to address this question. 

5. Summary 

Two fermionic RG methods were developed over the last years to investigate the physics 
of inhomogeneous LLs, with a special emphasis put on the transport properties of such 
systems. The poor man's RG (resummation of "leading-log" divergences) and the fRG 
(full functional dependence of the effective impurity potential) were set up to supplement 
the results obtained using bosonization for an effective field-theoretical model (the local 
sine-Gordon model) in various respects. We critically reviewed the results obtained 
within both approaches for the simplest setup, a ID quantum wire with two-particle 
interaction and a single local impurity. This system is characterized by two fixed points: 
the "perfect chain fixed point" which is unstable for repulsive two-particle interactions 
and the "open chain fixed point" which is stable. Close to the fixed points observables 
show power-law scaling as a function of energy variables. In the literature both fermionic 
RG methods are described as being applicable to capture the full crossover from one 
to the other fixed point of the single impurity problem while keeping interaction effects 
only to leading order. 

By comparing the outcomes of both approaches and relating them to what is known 
from bosonization we collected evidence that the fRG can be used to describe the 
crossover while the poor man's RG seems to have difficulties close to the "perfect chain 
fixed point" . Comparing the generalized Breit-Wigner forms in which the conductance 
in both approaches can be written we argued that in the poor man's RG the derived 
fiow equation for the transmission ai kp must be supplemented by at least a second 
independent equation. For a single impurity this does not have obvious consequences 
for the scaling behavior of the linear response conductance in the limits of weak or strong 
bare impurities because of a subtle relation (which can be traced back to what is called 
"duality" in the bosonized model) between the scaling exponents of the two fixed points. 
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For weak bare impurities and on large to intermediate energy scales the renormalized 
impurity potential decays with an exponent K within the fRG (which was argued to be 
consistent with what is found in bosonization) while it seems to decay with the exponent 
1 within the poor man's RG. In the fRG approach one obtains a one-parameter scaling 
function which depends on the strength of the two-particle interaction while the poor 
man's RG leads to the noninteracting scaling function l/(l-|-a;^). These three differences 
can be traced back to the main difference in the construction of the two fermionic RG 
methods. In the poor man's RG only the kp —kp (and vice versa) scattering is 
kept, while the fRG includes the mutual feedback of all scattering channels. In fact, 
for a single impurity the poor man's RG equation (|22l) can be obtained from the fRG 
approach by the following steps: (i) One neglects the flow of the two-particle vertex, (ii) 
One considers scattering states instead of Wannier states, (iii) The plane wave matrix 
elements of the two-particle interaction are assumed to be constant (up to the g-ology 
classification), (iv) In the overlaps of the plane wave and scattering states one assumes 
the asymptotic form of the scattering states even close to the local impurity, that is the 
nontrivial form of the scattering states generated during the RG flow is neglected. This 
last step leads to the decoupling of the scattering channels. Details on this are given in 
the Appendix. 

Both approaches can be extended to study more complex situations such as 
junctions of several wires, to characterize the corresponding low-energy physics. In this 
respect the bosonization approach is less flexible. It is interesting to address the question 
how both approaches compare in situations with a richer fixed-point structure in which 
the scaling exponents are not related by a simple transformation (such as K ^ 1/K 
for the single impurity case). A more complex behavior is found for Y-junctions pierced 
by a magnetic flux as studied in Refs. [3TJ using field-theoretical methods, in Ref. [15] 
applying the poor man's RG, and in Ref. [24j by fRG. In the field-theoretical approach 
it was possible to identify one of the several fixed points found using the fermionic RG 
methods and determine its scaling dimension as a function of the LL parameter K. This 
fixed point is characterized by a maximal breaking of time-reversal symmetry, that is 
the transmission probability from leg i to leg j (with i,j = 1,2,3) of the Y-junction 
is unity if i > j and zero if j > i (or vice versa). For this fixed point the scaling 
dimension of both fermionic RG methods agree to leading order with the one found for 
the field-theoretical model. This has to be contrasted to the situation at another fixed 
point, which was so far not identified in the field-theoretical approach. Both fermionic 
approaches predict that at this fixed point the transmission probability between all legs 
is equal to 4/9 ("perfect junction fixed point"), while the scaling exponents obtained 
differ even to leading order in the two-particle interaction. This situation might present 
an example in which the discussed problems of the poor man's RG have a more severe 
effect than in the case of a single impurity. To obtain more insights on this it would 
be very desirable to identify and characterize the "perfect junction fixed point" of a 
Y-junction using a method which does not require approximations in the strength of 
the two-particle interaction. 
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In a very recent preprint Aristov and Wolfle [32] derive a RG flow equation for 
tfie conductance of a LL with a single impurity by keeping more than the "leading-log" 
contributions. To describe the underlying homogeneous LL the TL model is used (as in 
bosonization). This approach seems to overcome the problems of the poor man's RG 
and in addition can be used for arbitrary amplitudes of the TL two-particle interaction. 
It would be very interesting to see if this method can be extended to treat systems with 
spin and with more complex inhomogeneities. 

Note added in proof. After submission of the present manuscript we were 
informed by D. Sen that he redid the poor man's RG calculation [15j for the Y- 
junction. He now obtains a scaling exponent of the "perfect junction fixed point" which 
is consistent with the result derived using fRG [23]. 
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Appendix A. Deriving the poor man's RG equation from the fRG 

The poor man's RG equation ( l22i) can be derived within the fRG framework. To achieve 
this we slightly modify the steps which led to Eq. ([7]). For the interacting wire we 
consider the spinless ID electron gas model instead of a lattice model. Furthermore we 
do not integrate out the "adiabatically" connected noninteracting (ID electron gas) leads 
before setting up the fRG procedure. As a consequence we have to deal with infinite 
instead of finite matrices. If the approximations which led to Eq. ([71) are applied this 
does not change the formal structure of this flow equation for the frequency independent 
self-energy. We now neglect the flow of the two-particle vertex altogether and replace 
F"^ by the antisymmetrized interaction v. The only energy scale which appears in the 
poor man's RG is the cutoff A. In the fRG we therefore consider the T = limit and 
later introduce temperature by stopping the RG flow at A = T. After these steps the 
approximate flow equation for the self-energy in the momentum state (plane wave) basis 
reads 



with the plane wave matrix elements Vk',q';k,q of the antisymmetrized two-particle 
interaction. Within the static approximation Q^{iu) can be interpreted as the resolvent 
[iu — Hq — Ti^]~^ of an effective single-particle problem characterized by the single- 
particle version of the noninteracting Hamiltonian ho and the flowing scattering potential 
S^. We emphasize that the effective potential is restricted to the finite interacting 
part of the system. At fixed A we can therefore define scattering states \k,+)j^, with 




(A.l) 
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\k, +)^ = Ek |^+)a = hQ + T,^. For large |x| the scattering states are given by 

{k > 0) 



gjfcx _|_ ^Ag ikx ^ 



2n 



— OO 

for X ^ OO . 



(A.2) 



A corresponding expression holds for < 0. The scattering states can be expressed in 
terms of the plane wave states using the T-matrix [33] 

1 



1 + 



\k) . 



(A.3) 



Ek- ho + iO' 

The transmission and reflection are given as plane wave matrix elements of the 
T-matrix 



i^k 



l-—Ak\T\eu + iQ)\k) 
2Tii 



(-A;|T^(£, + zO)|A;) 



\Vk\ 



(A.4) 



with the velocity = de^/dk. Using the distorted wave Born approximation [33] one 
can show that a change of the scattering potential + dT,^ leads to the changes 



\Vk\ 

2Tri 

\Vk\ 



'A ' 



,{-k,-\d^^\k,+)j 



of the transmission and reflection amplitudes. With the relation 



I^,-)a =tV\k,+)^ + r^:\-k,+) 
between the scattering states one obtains 



(A.5) 



(A.6) 



dti 



2m 

\Vk\ 
2TTi 



{tt A{k,+\ di:^ \k, +)^ + r^, A{-k, +1 rfS^ \k, +) J 



drt= ---{r^ A{k,+\d^''\k,+)A + t'^kA{-k,+\d^''\k,+)^) . (A.7) 
l^fc I 

We note in passing that these expressions also hold for lattice models. Combining 
Eqs. f lA.ip and flA.7p and transforming into the scattering state basis gives 



dtt 
dA 



drt 



I 

\vk\ 



dA \vk\ 



^aJ 

E [dq-^(rt 



A -s,A I A -s,A 
k '^k,q;k,q ' '^-fc '^-k,q;k,q I ' 



A -s,A I ,A -s,A 

'^k,q;k,q ' ^ -k '^-k,q;k,q I ' 



with the A-dependent scattering wave basis matrix element 



^s,A 
'^k'V±I 



j dqdq'dpdp' {k', + \q') {I', + \p') {q \k, +) {p \l, +) (A.9) 

of the antisymmetrized two-particle interaction. In the following we focus on the first 
equation in (lA.Sp . 
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To compute "W^;^/.;, i we implement further approximations. In a first step we project 
the plane wave matrix element Vq'^p'-q^p on the two Fermi points ±kp, that is we consider 
the g-ology model [29] with 



- VpOi S{q' + p' — q — p) 

{1 if sign q' = sign q = — sign p' = — sign p 
—1 if sign q' = sign p = — sign p' = — sign q (A.IO) 
otherwise 

and a as defined in Sec. [U The intra-branch interaction (usually called in the g-ology 
notation) only leads to a renormalization of the Fermi velocity. With respect to the 
exponent a of the power-law scaling this is an effect of higher than linear order in the 
interaction and can thus be neglected. Secondly and most severly, in computing the 
overlaps of plane wave and scattering states we take the asymptotic form Eq. (\A.2\i of 
the latter for all a; ^ 0. We thus neglect the nontrivial structure of the scattering states 
within the interacting region. The largest contributions to the integrals in Eq. flA.Op 
follow from momenta close to ^kp. In {k,+ \q) terms with denominators k ± q + iO 
appear. Depending on the signs of k and g as a third approximation we neglect those 
terms in which the sum or difference of the momenta becomes of order kp while the 
terms with small denominators are kept. For k,q > this gives 

j_ ( 1 4* 



{K+\-q) 

{-K+\q) 



2tt -k + q + iO 

_^ ' -k 

27r -A; + g + iO 



[-k,+ \-q) ^ — {- + , . (A.ll) 

' ' ' 2n\k-q + i0 -k + q + iOj ^ ^ 

As a last approximation we linearize the dispersion Sk = ±VF{k =F kp) appearing in 
Eq. (jAilD. 

Under these assumptions a lengthy but straightforward calculation yields 

# = tti^ - \tt\') • (A.12) 

dA Vk ^ ^ A-iVFik-kp) 

Remarkably, the different scattering channels are now decoupled, that is in computing 

the transmission at momentum k only the same k enters on the right hand side of the 

flow equation. At the Fermi momentum one obtains 

or 

dt^ 

din(A;i;./^) -<'^-i'tr) (a-m) 

which is the same as the poor man's RG equation 
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